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Abstract 

We derive general rigorous lower bounds for the average ground state 
energy per site e''*' of the quantum and classical Edwards- Anderson spin- 
glass model in dimensions d — 2 and d = 3 in the thermodynamic limit. 
For the classical model they imply that e'^' > -3/2 and e*^' > -2.204 • • •. 

1 Introduction 

The study of ground states (GS) of short-range spin-glasses is a very difficult 
subject, which is demonstrated by the fact that finding the GS of a spin-glass 
in a three-dimensional lattice is an NP-complete problem [T: these systems 
are thereby placed among those with the highest computational complexity. 
After several exact results for finite lattices have been reported (see later), a 
fundamental problem poses itself : how far removed are such results from the 
thermodynamic limit? We propose to offer a first (partial) answer in this paper. 

Our idea is to apply a method of Bader and Schilling [2] (for short BS 
method), which yields exact lower bounds to the GS energy per site for quan- 
tum Hamiltonians in the thermodynamic limit in terms of known results for 
finite lattices (cells), to random systems, such as the quantum anisotropic ver- 
sion of the Edwards- Anderson (EA) [3] spin-glass model treated in [HI, which 
includes the EA model as a special case. Our main observation is that the BS 
method is very naturally adapted to the concept of independence (of the cou- 
plings) in contrast to the standard methods of obtaining lower bounds, such as 
Temple's inequality 0], where an expectation value of the Hamiltonian appears 
in a denominator. Note that, in some instances, lower bounds given by the 
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BS method may be exact: examples occur in the theory of quantum spin sys- 
tems with competing interactions ([2],[5]), as well as in a theory of anisotropic 
ferromagnetic quantum domains [5], but in the present case this is not to be 
expected due to the aforementioned complexity. For two nice textbook accounts 
of rigorous methods in the theory of disordered systems, see [53] and pS] , 

Most of the numerical research on the EA spin glass ended at about 1980: 
we refer to Binder's review article [7] for the results of Monte-Carlo simulations, 
which predict a GS energy per lattice site of about —1.4 for the two-dimensional 
model, and about —1.9 for the three-dimensional model, but, as far as we know, 
no exact results exist, except for finite lattices, which we shall refer to shortly. 
Interest remained, however, in connection to the possible applicability of Parisi's 
theory to the EA spin glass : see [15] for a recent reference. Concerning lower 
bounds, Hoever and Zittarz |25j report a surprising value of —1.4219, in a re- 
markable paper which, we believe, merits to be revisited. Although they refer 
to this bound as "rigorous" , we believe that it is not, at this point, for a number 
of reasons: it would be so if the conjectured analyticity of their expansion held 
up to T = 0, but this is an open problem. Otherwise, it seems that a more 
careful estimation of the rest term after a given finite number of terms would 
be necessary, as well as a proof of the interchangeability of the limits T — 
and the thermodynamic limit in this expansion to finite order; this problem also 
occurs in connection with the ground state (or residual) entropy [TS] . There are 
also important bounds due to Kirkpatrick [22| . which will be commented upon 
in the conclusion. 



2 Results 

As in [5], for each finite set of points A C Z'^, where the dimension d will be 
restricted to the values d = 2 and d = 3, consider the Hamiltonian 

^fA({J}) = ^ J.,,<i>.,, (l-a) 

where with i,j G A are self-adjoint elements of the algebra generated by 
the set of spin operators, the Pauh matrices trf , af , trf , « G A, on the Hilbert 
space Ha = ^ieACf, given by 

^^,3 = a^a^a^ + ayafa^ + a.a^a] (1 - b) 

for |i — j| — 1, and zero otherwise. The random couplings J^j, with |i — j| = 1, 
are random variables (r.v.) on a probability space (fi, J^, P), where is a sigma 
algebra on n and P is a probability measure on We may take without 

loss of generality 

n=XBdS (2 -a) 

where 5 is a Borel subset of R, is the set of bonds in d dimensions, and 
assume that the Ji j are independent, identically distributed r.v.. In this case, 
P is the product measure 

dP = Xg^dPo (2 - b) 
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of the common distribution Pq of the random variables, which will be denoted 
collectively by J. The corresponding expectation (integral with respect to P) 
will be denoted by the symbol Av. We have to assume that 

Aw(J,j) = (2-c) 

for all «, j G Z"*, i.e., that the couplings are centered. Let Ea denote the GS 
energy of i?A, i.e., i?A = inf spec{H\). The following result was proved, among 
several others, in [9]: 

Theorem 1 (^) For P- a.e. {J}, the limit below exists and is independent 
of the b.c: 

^ ^ (3 - a) 

Ay'Z" |A| 

and 

where |A| denotes the number of sites in A. Finally, 

e'-"^ > 6^"+'^ (3-c) 

(3-a) is the far-reaching property of self-averaging (see [TO] for a discussion) 
from which it follows that, P- a.e., 

e(^) = linr (3 - d) 

AyZ'i |A| 



Let An denote a square with N sites if rf = 2 or a cube with N sites if d = 3, 
and write H^\j) = Ha,„{J)- We now adopt periodic b.c. for simplicity, but 
the final result is independent of the b.c. due to theorem 1. We may write 



where is given by 

Above, Cd are factors which eliminate the multiple counting of bonds, i.e., 

C2 = 1/2 (4 - c) 

and 

C3 = 1/4 (4 - d) 

(2) 

and A„ is a square labelled by a site n, for which we adopt the convention, 
using a right-handed {x,y) coordinate system, that n ~ {nx,ny) is the vertex 

(3) 

in the square with the smallest values of Ux and Uy . Similarly, An is a cube 
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labelled by a site n = {nx,ny,nz) with, by the same convention, the smallest 
values of n^, % and riz- Due to the periodic b.c, the sum in (4-b) contains 
precisely N lattice sites. The notation is short-hand for its tensor product 
with the identity at the complementary lattice sites in A;v\An'^^ Let Eq^''^\j) 
denote the GS energy of H^n\J) , e''^''^\j) the GS energy of Hi^\j) and 

4'^)^M4"''^)(J)) (4-e) 

By the condition of identical distribution of the r.v. Ji^j , e'q^ does not depend 

on n, which is implicit in the notation used. We have 
Theorem 2 The following lower bound holds: 



Proof. By the variational principle 

inf 



4^''^)(J)=,,inf {^,h(^\j)^1,)> 



J2 , inf (^P,Hi^\j)^) 

■-r (V'>i/')=i 



> 

(5-6) 

Let (V'a) denote the normalized eigenstates of H^\j) corresponding to eigen- 
values {eL"'^'>} with < a < and ^"'''V) < • " • eI"'^\J) J™ + 1 = 2l^"'*l 
denoting the dimension of the Hilbert space over We may write the 

general normalized wave- function tp in = 'H^w ^ ^AnVa''*' form 

tp = J2a p Ca,fiipa <8> with J^a l^ct,^]'^ = Ij and {(j)^) an orthonormal basis 
of ^AnVa'"*'" ^^^^^ H^\j) equals the identity on the latter sites, 

inf {^,Hi^\m = 
(V,V')=i 



inf^|c„,^p£;("''^)(J) 
>E^''''^\j) 



(5-c) 

which, together with (5-b), yields 

neAiv 

We now remark that , when performing the average over the J on the r.h.s. of 
(5-d), each r.v. Eq"''^\j) is the product of a function of the couplings in An^ 
and the identity function of the couplings in Ajv\a1''^; by mutual independence 
of the Jij this average equals, for each n, an average restricted to the couplings 
in A^'^'' which equals the left hand side of (4-e): this proves (5-a). q.e.d. 
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3 Application to the classical Edwards- Anderson 
spin glass 

The classical EA spin glass is the special case ax = ay = in (1-b); we also 
set ttz = 1. In this case e^'^\ given by the r.h.s. of (3-d), is invariant under the 
"local gauge transformation" erf —erf together with Jij —Ji.j'^j\\j — 'i\ = Ij 
whatever the lattice site i. Accordingly [3T] an elementary square (" plaquette" ) 
P is said to be frustrated (resp. non-frustrated) if Gp = Yl(^i j)^p Ji.j ~ —1 
resp. -1-1. Note that Gp is gauge-invariant and that, for the quantum XY 
(or XZ) model defined by setting aj, = in (1-b), e^'*-' is also locally gauge- 
invariant if we add to the above definition the transformation erf — )■ —erf, i.e., 
the transformation on the spins is defined to be a rotation of tt around the y - 
axis in spin space. 

Since the Pauli z-matrices commute, finding the minimum eigenvalue of (1- 
a) in the classical EA case is equivalent to find the configuration of Ising spins 
ai = ±1, denoted collectively by cr, which minimizes the functional 
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The minimal energy of a frustrated plaquette equals Epj — —2 and of a non- 
frustrated plaquette £'p.„ — —4. We have, now: 

Proposition 1 Let dPo in (2-b) be the Bernoulli distribution dPo — l/2{Sj+ 
6-j) and set for simplicity J — 1. Then: 

e(2) > _3/2 (6 -a) 

and 

,o^ , ,36096 
e(3) > _i/4 6-6 

- ' 4096 ^ ' 

Proof. For d — 2, there are exactly eight frustrated plaquettes (with Epj- = 
—2) and eight non-frustrated plaquettes (with Ep^u = —4), hence by (4-b), (4-c) 
and (4-e), E^^^^ = (1/2) (1/16) (-2 x 8 - 4 x 8) = -3/2 and (6-a) follows from 
(5-a). For d ~ 3, there are 2^^ possible configurations of the J and 2* possible 
values of the products CT^crj, hence a total of 921600 parameter values, almost 
one million, indeed huge in comparison with 256 for d — 2. We obtain the r.h.s. 
of (6-b) by (4-b),(4-d) and the use of a convenient program which computes 
£'q^-'/c3, with Eq^'' given by (4-e). In order to do so, the program performs a 
finite number of operations involving just sums, subtractions and comparison of 
integers, hence the final result is exact, q.e.d. 

Remark 1 Due to (3-b), the approach to the thermodynamic limit is in a 
monotonically decreasing fashion. Thus any exact value for the ground-state 
energy per site for a finite sample is an upper bound to e^''^ and, by (3-c), any 
upper bound to e*^^-* is an upper bound to e*-^-*. It is, however, difficult to ascer- 
tain whether some of these "exact values" are rigorous. Those obtained by the 
so-called branch-and-bound algorithm are rigorous, because rigorous numerical 
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error bounds of 0(1/L), for a. L x L lattice, in two dimensions, may be derived: 
see "Exact Ground States of Two-Dimensional ±J Spin Glass", by C. De Si- 
mone, M. Diehl, M. Jiinger, P. Mutzel, G. Reinelt, and G. Rinaldi, available 
at the web, in particular their Proposition 1. A similar rigorous result |19j (we 
thank Prof. S. Homer for this reference) is also available at the web, yields the 
GS energy per site for a 5X5X5 lattice. Hence, by monotonicity and according 
to the latter result, 

e(3) < -1.759 (7 -a) 



and in view of the former. 



e 



:2)<-1.39--- (7-5) 



where the number on the r.h.s. corresponds roughly to the best result estimated 
from their figure. An early reference to the case d = 3 is [11] . 

Remark 2 In [13], B. Derrida introduced the now famous random energy 
model (REM) as an approximation to the EA spin glass. In the case of the 
Bernoulli distribution, he obtains the lower bounds (see his eq. (54)): 

e^^^ > -1.560 (7-c) 

and 

e(3) > -1.956 {7-d) 

which are rigorous for the REM. It is interesting to remark that (7-c) compares 
well with (6-a) and (7-c) also compares well with (6-b), which may be written 

e(3) > „2.204--- (6 -61) 

Remark 3 To the best of our knowledge, proposition 1 provides the first 
(nontrivial) rigorous lower bound both for e^'^^ and e*^^^. Using the natural misfit 
parameter 

\E"^\ ^ ' 

as a measure of plaquettes frustrated or bonds unsatisfied (see (4) of [H]), 
where Eq denotes the ground state energy of the frustrated system and E"^^ is 
the ground state energy of a relevant unfrustrated reference system, we find 
from proposition 1 in the d = 2 case the lower bound m > 0.25 and for d = 3 
the lower bound m > 0.26 • • ■: thus, in both cases, the measure of frustrated 
plaquettes or unsatisfied bonds as defined above is at least of the order of 1/4. 



4 Conclusion 

Consider the dual lattice, formed by the centers of the squares {d = 2) or cubes 
{d — 3) of the original lattice. Toulouse ST considered paths (strings) made up 
of lines joining neighboring points in the dual lattice and crossing the unsatis- 
fied bonds, and suggested that the GS is characterized by the configuration(s) 
minimizing the average length L of such strings (in two dimensions, in three 
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a much greater complexity arises and the average surface-to-perimeter ratio of 
"minimal covering surfaces" 7 (see [22]) is the quantity to consider). Based 
on this suggestion, Kirkpatrick [3^, in a beautiful paper, obtained the bounds 
(6-a) and 

e(3) > _2.25 (9) 

under assumptions which seem to be eminently reasonable (L > 1 and 7 > 1/4) 
but are actually very hard to prove. Since then, several suggestions that the 
emergence of certain global configurations (i.e., configurations of the lattice as 
a whole), such as percolation of zero energy configurations (see [5], [10]) or 
"domain walls" associated to low-energy excitations consisting of clusters of 
overturned spins (with respect to some ground state) [33], with certain scaling 
characteristics, are responsible for the GS and low-temperature properties of 
spin glasses. We refer to [M] for a more complete list of references up to 1987. 
The fact that the r.h.s. of (6-bl) is even slightly better than (9), which relied 
on heuristics based on a global property, suggests that local properties might 
play an equally important role in the low-temperature physics of short-range 
spin glasses. It must be emphasized that for d = 3 there exist local correlations 
between the frustrated plaquettes, for instance, for d = 3, the product of the 
Gp over all plaquettes in the unit "cell" A^f^ equals +1, therefore there must be 
an even number of frustrated (and unfrustrated) plaquettes in a cell configura- 
tion. Therefore, unlike the case d = 2, not every plaquette configuration occurs 
independently in one cell configuration. 

We finish this paper with two remarks. In a low temperature expansion, 
the first term in the partition function is the GS degeneracy. Since Parisi's 
seminal work [T3], it is assumed that this degeneracy becomes infinite in the 
thermodynamic limit, and, indeed, it was recently shown numerically |15| that 
for the EA spin-glass with Bernoulli random couplings the overlap introduced 
by Parisi |14j is a good order parameter. In this connection, there exist rigorous 
lower bounds for the entropy density at zero temperature for d = 2 and d = 3: 
the first was due to [l7], which was greatly improved by [18]. For d = 3, see 

The second remark concerns the quantum model (1), to which the method 
is applicable: for ctj^ = we have the (quantum) XY model, which may define a 
different universality class, and it would therefore be very interesting to see how 
the present bounds would change with (fixing = 1). It is also, conceptually 
as well as practically, important to view the Ising model (1-c) as an anisotropic 
limit (e.g., as a^; and aj, — ;> in (1-b)) of quantum models: indeed, the 
reason why the critical exponents of the Ising model in three dimensions (see, 
e.g., [28]) are so close to those measured in real magnetic systems is that most of 
the latter are highly anisotropic. Thus we feel that it is important that general 
proofs are stable by small quantum perturbations. In addition, the only natural 
(i.e., not imposed) dynamical evolution is the quantum evolution. One exact 
example of the latter is furnished by the model treated in [1^ , to which we also 
refer for a discussion (with references) of the role of probability distributions in 
random systems. 
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